Dynamical properties of the pinned Wigner crystal 
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We study various dynamical properties of the weakly pinned Wigner crystal in a high magnetic 
field. Using a Gaussian variational method we can compute the full frequency and field dependence 
of the real and imaginary parts of the diagonal and Hall conductivities. The zero temperature Hall 
resistivity is independent of frequency and remains unaffected by disorder at its classical value. We 
show that, depending on the inherent length scales of the system, the pinning peak and the threshold 
electric field exhibit strikingly different magnetic field dependences. 
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Though the exotic possibility of electron crystallisation 
was discussed decades ago by Wigner ^ , its experimen- 
tal realization has been a challenge due to difhculties in 
obtaining sufficiently low density electron systems. How- 
ever, this problem can be circumvented by subjecting the 
2DEG to large magnetic fields which facilitate crystalli- 
sation of even dense electron systems. The quest for the 
Wigner crystal (WC) in mono ||-||] and bilayer quan- 
tum Hall samples indicated the existence of a quaint in- 
sulating state at filling fractions where crystallisation was 
theoretically expected. These and other detailed studies 
[0 of the insulating state revealed that the diagonal resis- 
tivity Pj-j. diverges as the temperature T — > and shows 
activated behaviour at finite T whereas the Hall resistiv- 
ity p^y is temperature independent and has a value close 
to the classical Hall value. 

Although measurements of activated linear and non- 
linear dc conductivity and the luminescence spectrum 
of radiative recombination were consistent with in- 
terpretations in terms of a pinned WC, the finite value 
seen for p^y was unexpected. This prompted other in- 
terpretations of the observed insulating phase P,p^ and 
in particular the existence of a new phase, the Hall insu- 
lator (HI) llfl]. The Hall insulator is defined as a phase 



where limt^^o 



cste which yields pxx oo and 



a finite pxy in the limit T ^ 0,lu ^ 0. This was proved 
only for non interacting electrons in a random poten- 
tial (i.e an Anderson insulator in presence of a magnetic 
field) and qualitative arguments suggested that it holds 
for interacting systems as well. However none of these ar- 
guments take into account the possible local crystalline 
order which could result in radically different physics as 
compared to the disordered electron fluid. Indeed period- 
icity plays an important role in other disordered systems, 
such as vortex lattices [ pl] |. 

It is thus of prime importance to investigate in de- 
tail the transport properties of a pinned WC. Transport 
properties are especially important here because of the 
extreme difficulty of a direct experimental verification of 
local crystalline order. One of the few theoretical at- 
tempts made to predict these properties was that of Ref. 
| p2[ where the related problem of charge density waves 
(CDW) in a magnetic field was studied. The harmonic 



approximation used, however, did not allow the extrac- 
tion of the detailed frequency dependence of the conduc- 
tivities. Later works focussed on the sliding state 
and the effects of free carriers, or on the effect of strong 
disorder [Q. A noteworthy point is that none of these 
calculations consider both a lattice structure and modu- 
lation of disorder at scales smaller than the lattice spac- 
ing . This feature which is absent in CDW turns out 
to play a crucial role in the physics of the WC. 

In this Letter, we compute for the first time the real 
and imaginary parts of the frequency dependent conduc- 
tivities of a weakly pinned WC. We find that even some of 
the main features derived in Ref. [Q are incorrect. The 
results we obtain provide a basis for comparison with 
recent experiments which map out the low frequency be- 
haviour of the conductivity [ p^p6| . 

Our starting point is the WC in a magnetic field B 
with lattice spacing a modelled by an elastic hamilto- 
nian ||l2|. The electrons at site i are displaced from 
their mean equilibrium positions by u(Ri,i). We 
also take into account the Coulomb repulsion between 
density fluctuations. We use the following decomposition 
u(q) — quL(q) + qAzuT(q) where L, T denote the longi- 
tudinal and transverse components. The corresponding 
action in the imaginary time formalism is 
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Pm, Pc are the mass and charge densities, c and d are 
the shear and bulk modulus respecitvely. For the WC, 
the presence of Coulomb forces results in a bulk 

2 

modulus d — — much greater than the shear modulus 

2 

c ~ ^ (eo is the dielectric constant of the substrate). 

— PcB I Pm is the cyclotron frequency and the Matsub- 
ara frequencies at temperature T are ujn = ^irn/ ph where 
/3 = 1/T. From (|), we see that the Coulomb interaction 
affects only the longitudinal modes. The magnetic field 
couples the transverse and longitudinal modes. ( ) 
denote averages over quantum and thermal fluctuations 
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and are disorder averages. For the pure system the 

quantum fluctuations result in (u^) where Ic is the 
magnetic length Ic — y^h/eB. The time-independent dis- 
order potential VF(x) is short range correlated (of range 
rj) W{x)W{x.') = A6rf{:xi — x') and couples to the den- 
sity of electrons p{x, t) — J2i — Ri — u(Ri, t)). Using 
the decomposition of the density into lattice harmonics 
(valid in the absence of topological defects) and repli- 
cas to average over disorder we obtain the effective action 
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J2 cos(K.(u''(xt) - u''(x, r'))) (2) 

a,b,K 

a, b denote the replica indices, K are the reciprocal lattice 
vectors and Ak ~ Aexp— K^r^. Note that it is impor- 
tant to retain all harmonics. The disorder averaging also 
yields a term quadratic in the displacements which has 
been absorbed by a shift in u. This shift does not affect 
the conductivity and we neglect it henceforth. Actions 
similar to (||) can be used to describe 3-d classical prob- 
lems such as vortex lattices with correlated disorder 
and long range interactions. 

To study this model we use the gaussian variational 
method (GVM). This quantitative method allows us to 
compute the Green functions and hence the conductiv- 
ity of the system described by (|^). Unlike previously 
used methods the GVM is self-consistent, has no un- 
determined adjustable parameters and also incorporates 
important physical feautures of the problem such as the 
existence of many metastable states. It allows one to go 
beyond simple static arguments as will be shown later. 
We introduce the variational action [flSl 



^0 



(3) 



where the Green functions G"^ are the variational pa- 
rameters ( a,/3 L^T and summations over repeated 
indices are implicit). They are determined by solving 
the self-consistent saddle point equations obtained by ex- 
tremizing the variational free energy i^var = + {ScR ~ 
50)0- The method extends the one used in ||l8| and all 
technical details will be presented in |l9j. As in the 
solution has a replica symmetry broken structure nec- 
essary to correctly describe the localization. The final 
result ||l9| is a closed set of equations for the connected 
part of the Green function G~^p — G~p which de- 
termine all physical quantities of interest here. These 
equations are respectively: 
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{cq^ +dq + PrnUJ^ + F) 



Gel = (c?' +dq + P„^ujI) + F + 
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{cq^ + p^ul + F) 



_ , (cg^ + P„,ojI + F){cq^ +dq + p„,wl -^ F) 

'^cLT ^ Pm^^n^c H 

Pm.l-^n^c 

with F — I{u)n) + S(l — 5n,o)- The localized phase is 
characterized by a non zero S, from which a length scale 
I can be defined through S = cl^^ . The function I{ujn) 
is defined as: 

2 /■'^ 

iM^-^j^ dT(l-cos(L^„T))(y'(B(r))-r(B)) (5) 

where the local diagonal correlation B{t) = 
l/2((w(0,r)-u(0,0))2) = 1/2{Bl{t) + Bt{t)) and the 
off diagonal part B are given by: 

Bt,l{t)^ ^ y X! GcT,l{^ - cos(w„t)) 



-^=75 / [2^iGcLiq,^n) + GMq.^n))] + 
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cq'^ + dq + 1] 
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Finally the set of equations close as S is itself determined 

by 
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(c<72 + dq + S)2 



(7) 



The primes denote derivatives. All information on the 
disorder is contained in the auxiliary function V[B] = 

In this paper, we focus on the transport properties but 
other quantities such as positional correlation functions 
can also be computed . The dynamical conductiv- 

ities are given by the standard analytical continuation of 
the Green's functions aapii^) = ip1^ujGap{q = 0,uj + ie). 
Rotational invariance combined with fact that a magnetic 
field breaks parity and time-reversal implies 
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iuj[—p,nUj'^ -I- S-f 
(S-p„C^2+J(^))2_^2^^2^2 

2 iujl-ipmUJUJc] 

(El - p„iw2 I{uj)Y - pI.uj'^ujI 
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In the absence of disorder one has / = S = in (||). 
Uxx vanishes in the dc limit = and has a 5— function 
peak at cyclotron frequency uj — lOc- On the other hand 
axy{ijJ = 0) = Pel B and a^y has a pole a,t uj = lOc- In the 
presence of disorder the crystal is pinned and conductivi- 
ties develop a new peak at the pinning frequency lo = ojp. 
Simultaneously, there is an upward shift of the cyclotron 
resonance peak from ojc by a quantity of order ujp. 
(4) To obtain the full frequency dependence of the conduc- 
tivities one needs to compute This can be obtained 
from the above equations which are valid for all values 
of B. Here we present a solution in the experimentally 
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relevant limit ^ dj c. A typical plot of Reaxx ob- 
tained by solving (j^) numerically is shown in Fig. 1. 
Since /(oj 0) = by definition in the pinned crystal, 
the dc value of cjxx is still zero but that of Gxy is zero in 
contrast to the pure case where it was finite. The peaks 
at the new resonance frequencies have a finite height and 
width due to disorder induced dissipation. The extent 
of this dissipation is determined by I{ujn) continued to 
real frequencies. Earlier results can be recovered by 
setting I{ijJn) = in all the equations. However, the 
presence of the I{LOn) term has many important physical 
consequences as will be discussed below. Firstly, in the 
absence of /(wn) the peaks would be delta functions at 
and ujc + with lu^ = Tj/uJc- In contrast, here the 
peaks are centered around a frequency ujp < and this 
shift is given by E. The peaks have a non-trivial struc- 
ture and are asymmetric about the resonance frequencies 
as can be inferred from (||) and seen in Fig. 1. This inval- 
idates the Lorentzian shape of the peaks which was used 
to arbitrarily broaden the delta functions in Ref. [|2| . We 
note that the peaks we obtain are much narrower than 
the Lorentizian broadened ones. 

For frequencies lo <^ iOp and ujp <ti lu ^ uJc, analytical 
solutions can be obtained. We find 



liuJn) = '\/2pmS H -==—\uJn\, 



UJ <t:. UJr, 



/K) = - log ^"^ " , 



(9) 



Using (||) in (||), we obtain the following low frequency 
behaviour for lo lOj, 



Reaxxi^^) = 
Imaxxii^) = P: 




Imaxyiio) ^ PcPrn 5 



In the region Up uj uJc '^G find using (^) 



RecTxyiuj) - 



(10) 
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Note that Reaxx and Reuxy are both quadratic in oj. 
Since the pinned WC has the characteristics proposed for 
the HI, it seems unnecessary here to invoke the existence 
of the HI as a new phase. The results of ( p^ ) can be used 
to calculate the dielectric constant e(w) = IvnUxxi^)/^- 
Its dc value is given by e = Pc/^- Thus the dielectric 
constant is also a measure of the characteristic frequency 
defined by disorder. 




FIG. 1. Reaxxii-o) (in units of hpm/ pc) as a function of to 
(in units of pc/Pm^) where h = cAa? /S.p'^ {e/Kf' for the case 
Tf < Ic for different values of the field B meausred in units 
of p^l pmh. The dashed line represents B = 0.33, the full line 



B = 0.4 and the long dashed line B 
magnification of the pinning peak. 



0.5. The inset is a 



Calculating the resistivities given by pafi = 
— o'q/3/(o'?2; + <^ty) find that the pinned crystal is in- 
deed insulating i.e., Pxx{^ = 0) = oo. More importantly, 
the Hall resisitivity pxy turns out to be independent of uj 
and T and has the same value as that in the pure system 
Pxyiiij) — B/pc- A similar result was argued to hold at 
T = in Ref. (l|]. At T > it would be necessary to 
go beyond the GVM approximation to ascertain whether 
Pxy still sticks to its classical value. Indeed the GVM 
misses soliton like excitations which are known to be im- 
portant for finite T physics |l8|j2^ ]. 

It is interesting to calculate the field dependences of the 
above quantities. This is of direct experimental relevance 
and we use our results to calculate the same. The field 
dependence of the pinning peak (whose width is naively 
of 0{ujp)) is governed by S whose value is in turn dictated 
by the relative sizes of the length scales rf and Ic- Tradi- 
tionally, in the context of CDW, S has been related to the 
Fukuyama-Lee length Ra at which relative displacements 
are of order a, as E ~ Ra^- However for the present case 
such a connection does not hold, because disorder can a 
priori vary at scales much smaller than a unlike in CDW. 
This is similar to the situation in vortex systems where 
pinning is controlled by the Larkin length i?c, defined as 
the scale below which the physics of can be described 
perturbatively by a model where uncorrelated gaussian 
random forces of strength Ay = act indepen- 

dently on each electron. Within this model Rc is given 
by ([u(R,) - u(0)]2)y = max[r2,l2] = When i?, » a 
the crystal is pinned collectively. In this regime using 
the GVM one finds S = cAR~^{a/^of. Here A is an 
overall constant, Ra pm.a'^/V^ and = max[rf,lc]- 
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This corresponds to E '--^ cR~^, implying that / ^ Rc 
(defined above), which shows that the length scale deter- 
mining the peak in the conductivity is Rc and not Ra- 
it is important to distinguish between these two lengths 
since Rc can have an explicit dependence on the mag- 
netic field. This yields two very different regimes. One 
is Tf < Ic, which gives S = bB^ {b = cAa^Ap~^(e/7i)^ ) 
hence ujp[B) oc and the pinning peak moves up and 
broadens with increasing field. This is the case illustrated 
in Fig. 1 where ReUxx has been plotted for various val- 
ues of B. The second regime is r/ > Ic leading to a E 
independent of B and LOp{B) oc B^^ . Thus the pinning 
peak moves towards the origin and gets narrower with in- 
creasing field. For CDW in a magnetic field ~ a and 
one is always in the second regime. In these two regimes 
the height of the pinning peak decreases as B~^ with in- 
creasing field. For the case rf ^ Ic, ojp does not increase 
indefinitively with B, rather there is a crossover to an- 
other regime when Rc ^ a where single particle pinning 
effects are dominant. Here the correspondence between 
Rc and I and hence S no longer holds. One finds then 

3 1 

that S cx i?2 and uip (x B^ provided Ic ^ Vf. In contrast 
the peak at uJc + always moves upwards with increas- 
ing B. A summary of the results in various regimes is 
given in Table ^ 

Another important measurable quantity is the thresh- 
old electric field Et necessary for the crystal to slide. 
This again shows the interplay between Ic and r f . Using 
collective pinning arguments |^ one gets Et = cRc^^o 
in the regime Rc S> a. The threshold field has the same 
regimes as above and the field dependence are shown in 
Table |. Note that for rf > Ic Et is independent of the 
field (as for CDW) . Since both Et and ujp are related to 
Rc one has Et oc Up but with a prefactor depending on 
and not on a as given by previous CDW estimates. When 
Rc < a, we enter the regime of single particle pinning. 
Using S oc i?2 ^ the threshold field is now Et oc TjIc oc B 
psi . Finally, due to the variation of S with the field the 
dielectric constant e will exhibit the behaviours shown 
in Table |. Therefore in addition to detailed frequency 
measurements of the conductivity, measurements of the 
B dependences of the dielectric constant could serve as 
an experimental signature for the WC. 

Some of the existing experimental results can be inter- 
preted within our theory. Contrary to previous estimates 
of the pinning frequency, it allows a scenario where the 
pinning frequency increases with the field as was seen 
in recent experiments JigI . A simultaneous increase in 
Et vs. B is observed which is in qualitative agreement 
with the above predictions. Our theory also predicts that 
the Hall resistance takes its classical value which is ob- 
served experimentally 0. However, many problems re- 
main both theoretically and in comparison with experi- 
ments. Experimentally the peak height in a{Ljj) seems to 
increase with B which we cannot account for at present. 
Some experiments m seem to report a different behaviour 



for the conductivity (see however |^J|). Finite tempera- 
ture effects also need to be understood since many exper- 
iments are performed in the d.c. limit at T > 0. Strong 
disorder effects have also to be understood. Both prob- 
lems require a careful treatment of the topological defects 
and solitons which are beyond the scope of the present 
study. While a phase transition similar to the one oc- 
curing in 3d vortex lattices |Tl| is unlikely in d = 2, one 
expects a marked crossover between a weakly pinned WC 
and a strongly pinned one. 

To conclude, we have developed a comprehensive the- 
ory for the WC pinned by weak disorder. In addition 
to detailed frequency dependences of the real and imag- 
inary parts of the conduct ivites, we have obtained the 
magnetic field dependences of various dynamical quanti- 
ties. We find that the magnetic field not only confines the 
electrons but also plays a crucial role in determining the 
response of the system to disorder. This dynamical effect, 
not captured by previous static approximations, allows 
the possibility of observing novel field dependences. 

We thank F.I.B. Williams for enlightening remarks. 
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Regime 


E 




Et 




r/ > Ic 




B-' 


B" 




rj < Ic 




B^ 


Bi 












Rc <a 




Bi 


B 





TABLE I. Magnetic field dependences of various dynami- 
cal quantities. 
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